

#taken from chanegpoint package R manual
# as used by Adam Kashlak in his significance paper February 2019 on academy award speech analysis
# (paper pdf in Kashlak reference on statswiki changepoint webpage)
#install.packages("changepoint")
set.seed(1)
x=c(rnorm(50,0,1),rnorm(50,5,3),rnorm(50,10,1),rnorm(50,3,10))
cpt.meanvar(x,penalty="Manual",pen.value="4*log(n)",method="BinSeg",Q=5,class=F)
plot(x)